Granular segregation as a critical phenomenon 
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We present the results of an experimental study of patterned segregation in a horizontally shaken 
shallow layer of a binary mixture of dry particles. An order parameter for the segregated structures 
is defined and the effect of the variation of the combined filling fraction, C, of the mixture on the 
observed pattern formation is systematically studied. The surprising result is that there is a critical 
event associated with the onset of the pattern, at C = 0.647 =b 0.049, which has the characteristics 
of a second order phase transition, including critical slowing down. 

PACS numbers: Ol.SOPa, 45.70Mg, 64.75.+g, 68.35.Rh, 45.70Qj. 
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Segregation is a counter-intuitive phenomenon where 
an initially mixed state of dry granular particles separates 
into its constituent components under excitation [Q. In- 
triguingly, it does not always happen and the conditions 
for its occurrence are difficult to predict. Segregation is 
not only of fundamental interest but it is also of practi- 
cal importance with applications in areas ranging from 
industry]^ to geology ||]. In recent years, there has been 
an upsurge of interest in small scale laboratory studies 
where vibration || , avalanching in partially filled hor- 
izontal rotating drums 0, 0] and stratification in verti- 
cally poured mixtures Q have all provided interesting 
examples of pattern formation. 

Our focus is on quasi-2-dimensional horizontally driven 
layers of binary mixtures of particles ||, as this 
gives the practical advantage that any collective behav- 
ior is readily visualized and gravity is effectively elim- 
inated. Moreover, the material is in contact with the 
drive throughout the motion. For this class of binary 
granular systems a qualitative segregation mechanism 
has been suggested [0 borrowing the idea of excluded 
volume depletion from colloidal systems and binary alloys 
[|l2|. In a driven 2D binary system of particles of differ- 
ent sizes, the packing fraction of the system is decreased 
if self-organized clustering of the larger particles occurs. 
Thereby considerably increasing the number of configu- 
rational states of the system. Hence ordered arrangement 
of the large spheres can increase the total entropy of the 
system by increasing the mobility of the small particles. 
This process has been referred to as entropic ordering . 
We believe that analogous attractive depletion forces are 
responsible for segregation in our granular system. We 
show that segregation undergoes a phase transition and 
occurs only for filling fractions above a critical value with 
the characteristics of a second order phase transition with 
critical slowing down. 

A schematic diagram of the top view of the apparatus 
is presented in Fig. |l|. It consisted of a horizontal smooth 
rectangular tray, of dimensions (x^y) = 180 x 90mm with 



a flatness of less than zb5/im, on which a binary mixture 
of particles was vibrated longitudinally. The tray, made 
out of aluminum tool plate, was mounted on a system 
of four high precision linear bearings and connected to a 
Ling electro- mechanical shaker with a servo feedback con- 
trol. The dynamical displacement and acceleration were 
monitored by a linearly variable differential transducer 
(LVDT) and two orthogonal piezoelectric accelerometers. 
The motion was checked to be unidirectional and sinu- 
soidal to better than 0.1%. The forcing parameters were 
kept constant at amplitude A = ±(1.74 ± 0.01 )mm and 
frequency of oscillation co = 12Hz. Static charging ef- 
fects were eliminated by coating the container's surface 
with a layer of colloidal graphite therefore making it con- 
ducting. Furthermore, provision was made to evacuate 
the apparatus and this had no influence on the segregated 
structures. 

The binary mixture consisted of poppy seeds {pps = 
0.20^cm~^), which were approximately ellipsoidal with 
major and minor axes d^^ = 1.07mm and d^J ~ dp~ ~ 



accelerometers 




FIG. 1: Top diagrammatic view of the apparatus. The exper- 
iment is set horizontally with gravity pointing into the page. 
The inset photograph shows an example of the homogeneously 
mixed initial conditions of the granular layer with phosphor- 
bronze spheres (White regions = 1) and poppy seeds (Black 
regions = 0). 
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FIG. 2: Series of space-time diagrams for increasing C. Each diagram is a stack, in time, of the mi-layer cut hne (LT in Fig. 
The LH side of the diagram corresponds to gaseous binary state. For larger values of C, on the RH side, segregated structures 
form within the first minute of forcing. All runs were started from homogeneously mixed initial conditions. 



0.50mm, respectively, and precision phosphor-bronze 
spheres {pps = 2.18gcm~^) with diameter dp^ = 1.50mm. 
The amounts used were always such that a shallow layer 
regime was maintained, with layer height h < 1.5mm. 

Forcing of individual particles arose principally from 
the frictional interaction with the base of the container 
and the effects of boundary walls were observed to be lo- 
calized. Although the motion of the container was sinu- 
soidal, the stick- and-slip nature of tray /particle contact 
induced stochastic thermalization; quasi-2-dimensional 
trajectories were obtained in x and with the character- 
istics of a 2D random walk [|l^ . The surface of the tray 
was accurately levelled with fine-pitch adjustable screws 
to within ±0.0086° and checked using a high precision en- 
gineers level. This ensured homogeneous forcing across 
the whole layer. 

The granular layer was imaged using a CCD camera 
mounted directly above the apparatus and 600 x 300 pix- 
els frames were digitized into a PC. Digital image pro- 
cessing techniques were employed to transform the ac- 
quired video frames into binary format, giving a value of 
(black) to pixels corresponding to regions of poppy seeds 
and a value of 1 (white) for those of phosphor-bronze 
spheres. Frames were thereby analyzed as 600 x 300 dig- 
ital matrices, allowing measurements on the dynamics of 
the granular patterns to be performed. 

We define a measure of the combined filling fraction of 
the layer, which we call the layer compacity, to be, 

xy 

where Nps and Np^ are the numbers of poppy seeds 
and phosphor-bronze spheres in the monolayer, Apg = 
0.87 ± O.OQmm^ and Api, = 1.77 ± 0.06mm^ are the two 



dimensional projected areas of the respective individual 
particles and x and y are the longitudinal and transverse 
dimensions of the tray. We regard C as the principle pa- 
rameter of the system. 

A digitized image of a mixture of the two types of par- 
ticles has been superposed on the schematic diagram of 
the apparatus in Fig. |l|. It is representative of the initial 
conditions which were consistently set using the following 
method. Firstly, Nps poppy seeds were vibrated at large 
amplitudes, A ^ ±bmm, creating an homogeneous and 
isotropic sub- monolayer. The phosphor-bronze spheres 
were then suspended, at a height h = 15mm above the 
poppy seeds layer, on a horizontal frame with a perfo- 
rated array of (57 x 28) holes, of 2mm diameter, arranged 
in a triangular lattice with a cell spacing of a = 2.75mm. 
A shutter was then opened and the 1596 phosphor-bronze 
spheres fell onto the layer of poppy seeds, creating a near 
homogeneous mixture of the two types of particles. 

The focus of the investigation was to systematically 
increase C, consequently decreasing the mobility of indi- 
vidual particles, and explore the conditions under which 
granular segregation took place. To achieve this, we in- 
crementally increased the number of poppy seeds, Nps, 
in the layer, in measured steps, with the amount of 
phosphor-bronze spheres, A/p5, held constant. We also 
experimented with changing the number of phosphor- 
bronze spheres and found that this did not qualitatively 
affect the main findings and the quantitative changes 
were found to be linearly scalable. We chose to perform 
the experiments by changing the numbers of poppy seeds 
since this provided finer steps in the control parameter, 
i.e. C = C{Nps). 

In Fig. H we present a series of space-time diagrams, 
for increasing C, which were constructed by sampling the 



3 




5-10 15 
Stripe width L (mm) 



CO 5 
0) 

2 4 





(C=0.873) 








* I^*(C=0.729) ; 










I^(C=0.606) 











50 100 150 200 
Time (sec) 

FIG. 3: a) Histogram of the L distribution (longitudinal 
stripe widths) for a single 600 x 300 pixel frame {C = 0.966). 
Solid curve is Gaussian fit. Dashed line is average value = 
(L) which we define to be the order parameter, b) Time series 
for the order parameter for various values of C A saturation 
level is attained after a fast initial segregation period. Solid 
curves are fits to Eqn. (0). 



mid-frame cut line along the granular layer's x-dimension 
(shown in Fig. [l] as the white longitudinal line LT), 
and progressively stacking, in time, 200 of these (600 x 
1 pixels) lines (acquisition rate of l/6Hz), for runs of 
lOmin. This period was found to be long enough for 
steady state conditions to be achieved. 

In the left-hand image of Fig. |, with C = 0.503±0.035, 
it can be seen that the binary layer remains mixed for the 
duration of the experiment. The system acts as a low 
density granular gas with two components and individ- 
ual particles describe quasi-2-dimensional Brownian-like 
trajectories. Similar behavior has been found in numeri- 
cal studies of two-dimensional binary granular gases with 
uniform thermalization in which a mixed non-equilibrium 
steady state is achieved at low filling fractions Q. The 
right-hand image in Fig. |, for C = 0.955 ± 0.081, dis- 
plays qualitatively different behavior. It can be seen that 
within the first minute of the start of the motion, small 
localized clusters of phosphor-bronze spheres (white re- 
gions) form and progressively coalesce with neighbouring 
ones. Thereby coarsening occurs, so that well defined 
stripes are eventually formed, aligned perpendicularly to 
the direction of forcing [^. 

At intermediate values of C, partially segregated states 
emerge and we aim to show that there is a critical de- 



pendence of the process on the layer compacity rather 
than a smooth emergence of the segregated structures. 
Thus there is a critical value, Cc, below which the binary 
granular layer remains mixed. Above Cc, as the layer is 
incrementally compacted, segregation structures develop 
in a nonlinear manner. Specifically, we treat the process 
as if it were a second order phase transition. 

In order to proceed, we first establish an order param- 
eter, ^, which we choose to be the characteristic lon- 
gitudinal width of the phosphor-bronze (white) stripes. 
For each acquired frame this was calculated by scanning 
through each of the 300 horizontal lines; the length of 
the white steps was measured, L(5zlJ), which typically 
yielded a distribution of ~ 6000 identified stripe widths. 
These were distributed normally, as shown in Fig. ||a, 
where a fitted Gaussian function has been superimposed 
on the experimental histogram of L. Since the distribu- 
tion is well defined, we have chosen the order parameter 
to be (/) = (L), i.e. the average value for the distribution 
of longitudinal widths of the segregated stripes. 

In Fig. we present time series of (j) for various C, 
where each run has been started from mixed initial con- 
ditions. At large compacities, e.g. C = 0.914 ± 0.077, 
there is a fast initial growth, eventually saturating af- 
ter ~ 1mm from the start of the experiment. At this 
point an approximately constant level is attained. The 
superimposed solid lines are fits of the form. 



a. exp 



t 



ts{C) 



(2) 



to the experimental time series of (j) where ^^(C), the 
value at which (j) saturates, is the segregation levels 
^(f){C) — a is the initial value, (j){t = 0), and ts(C), the 
time scale associated with the saturation of the order 
parameter, is the segregation time . 
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FIG. 4: Phase transition plot for segregation level. Solid line 
is theoretical curve for square-root singularity. Critical point 
occurs at Cc — 0.647 =b 0.049. The mixed phase exists below 
this point and the segregated phase above it. 
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FIG. 5: Plot of segregation time, ts, which diverges as the 
critical point Cc is approached from above. Solid line is theo- 
retical curve of the form: ts ^ {C — Cc)~^ ■ 



The C-dependence of the segregation level is presented 
in Fig. ^. We see that for low compacities, up to a 
value, Cc = 0.647 ± 0.049, the segregation level remains 
constant at (S^) = 3.21 ± 0.08mm. This corresponds to 
approximately two sphere diameters so that any domains 
are not classified as structures and the system is deemed 
to be in a mixed or unsegregated state. As C is increased 
past Cc, segregated stripes of increasing (j) emerge and the 
segregation level exhibits a square-root singularity. The 
solid curve in Fig. ^ is the line of best fit of the form. 



(3) 



where A is the transition scaling factor, {Cc^Ta(^(Cc)\ is 
the critical point at which the transition occurs and f5 = 
I is the order parameter exponent. The numerical values 
for the fitted parameters were A = 2.81 ± 0.07mm, Cc = 
0.647 ± 0.049 and ^^(Cc) = 3.21 ± 0.08mm. 

It is interesting to note that the generally accepted 
value of maximum packing fraction for a random closed 
packing for a 2D system of particles |l^ is around 0.82 
which is considerably larger than Cc for segregation in 
our system. This seems to indicate that the segregation 
transition occurs while the granular mixture is still in the 
monolayer regime. 

A further measure that supports the notion of critical 
behavior is the C-dependence of the segregation time, t^, 
extracted from the fittings of 0(C,t) to the function in 
Eqn. (^). Note that is is only defined for C > Cc, corre- 
sponding the region where segregation occurs. The mea- 
sured segregation time is plotted in Fig. |^ as a function 
C. As expected for a second order transition, is diverges, 
as C is decreased from above, near Cc This is usually 
referred to as critical slowing down where the divergence 
has the form of - (C - Cc)"^ [|l§. The solid line in 
Fig. H corresponds to the mean field approximation re- 
sult, 7 = 1 111. 



In summary, we have shown that granular segregation 
in a horizontally shaken layer displays both qualitative 
and quantitative features which are consistent with a sec- 
ond order phase transition, i.e. there is a critical layer 
compacity value below which mixing of the binary layer 
is guaranteed and above which reproducible segregation 
patterns form. The phenomena were not affected when 
the granular layer was evacuated, establishing that this 
is a fundamental process intrinsic to the dynamics of the 
granular layer. The critical exponents for the order pa- 
rameter and the segregation time agree with those of 
mean field type behavior, = ^ and 7 = 1, respectively. 
This is surprising since the system is far from equilib- 
rium; it is both forced and the principal energy losses are 
due to inter-particle inelastic collisions. A reconciliation 
between this results on granular segregation and recent 
theoretical reports of critical behavior in granular hydro- 
dynamics Jl^ and the existence of self-organized domains 
in granular compaction models jl^ may be crucial to our 
further understanding of segregation and other collective 
phenomena in granular materials. 
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